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We show that s-wave scattering resonances induced by dipolar interactions in a polar molecular gas 
have a universal large and positive effective range, which is very different from Feshbach resonances 
realized in cold atoms before, where the effective range is either negligible or negative. Such a 
difference has important consequence in many-body physics. At high temperature regime, a positive 
effective range gives rise to stronger repulsive interaction energy for positive scattering length, and 
weaker attractive interaction energy for negative scattering length. While at low-temperatures, we 
study polaron problem formed by single impurity molecule, and we find that the polaron binding 
energy increases at the BEC side and decreases at the BCS side. All these effects are in opposite to 
narrow Feshbach resonances where the effective range is negative. 



Polar molecular gas with strong dipolar interaction is 
a system of great interest to both atomic and molecular 
physics and condensed matter physics communities. Re- 
cent development of stimulated Raman adiabatic passage 
(STIRAP) technique has succeeded in creating a nearly 
degenerate gas of KRb molecules in its rovibrational 
ground state [T] . Unfortunately, the KRb molecules is not 
stable against two-body decay into K2 and Rb2 molecules 
[2]. Recently, theoretical calculation predicts that this 
problem can be avoided by choosing other combinations 
of alkali atoms, such as KCs, KNa, NaCs, NaRb and 
RbCs [3]. Applying STIRP technique to these species is 
now proceeded in many different laboratories all over the 
world. It is very promising that degenerate polar molec- 
ular gases can be realized within the next few years. 

For these molecules their electric dipole moments can 
be polarized by an external electric field (say, along z 
direction). As the applied electric field increases, the 
dipole moment d continuously increases from zero to the 
permanent dipole moment, and different molecules have 
different values of permanent dipole moment [3]. There- 
fore the two-body interaction potential studied here can 
be modeled as 

f rf 2 (1-3 cos 2 6) 

Vb(r) = ^ ' > rc (1) 

[00 , r < r c 

which contains a long range dipolar interaction and a 
short range non-universal potential. 9 is the angle be- 
tween f and z. Here for simplicity, we use a hard core po- 
tential to mimic short-range behavior of very complicated 
realistic potential, where r c is the core size. We can intro- 
duce a length scale D = md 2 /(2h 2 ), and then the typical 
dipole energy is given by h 2 D/(2m(r) 3 ) « (kpD)EY- By 
controlling the electric field, one can tune kpD from zero 
to the regime kpD s» 1 [3], i.e. the dipole energy is com- 
parable to the Fermi energy. 

Understanding the two-body problem of two strongly 
interacting quantum dipolcs is an important step toward 
understanding the rich many-body physics of this sys- 
tem. Previously, the dipolar interaction effects in the 
high partial wave channel have been studied extensively, 



for instance, the p-wave pairing [5J [B] and p-wave and 
d-wave Fermi surface distortion [5j [7] . While the dipolar 
interaction effects in the s-wave channel have not been 
well studied, because they always vanish in simple mean- 
field treatment. 

In fact, previous studies have shown that, as D in- 
creases, each partial wave channel will display a series 
of resonances where the scattering length of each par- 
tial wave channel diverges, and we call them dipolar in- 
teraction induced resonances (DIIR) [8HT0], Among all 
these channels, resonances in the s-wave channel are the 
most pronounced. However, both the two-body proper- 
ties and the many-body physics at DIIR have not been 
completely understood. In this letter we study both two- 
body properties of DIIR and two many-body problems of 
two-component fcrmionic molecules across a DIIR, that 
are, high-temperature regime with equal population, and 
single impurity in fully polarized Fermi gas. From these 
studies we want to emphasize the significant difference 
between s-wave DIIR and other s-wave resonances. Main 
conclusions are summarized as follows: 

(i) Near s-wave resonances the effective range ro is 
large and positive. This is qualitatively different from 
a magnetically tuned Feshbach resonance (MFR) stud- 
ied before, where r is either negligible (wide MFR) or 
large but negative (narrow MFR). Across different DIIRs, 
ro/D is a universal function of D/a s and r$/D ss 1.84 at 
resonance. Therefore &f?~o is of the order of unity and 
the effect of ro can not be ignored in many-body physics. 

(iii) At high temperature, the attractive interaction en- 
ergy of a DIIR is significantly smaller for negative scatter- 
ing length a s compared to a wide MFR; while for positive 
a s , in the upper branch, the repulsive energy of DIIR is 
much larger. Both can be understood by analyzing en- 
ergy dependence of the scattering length or the two-body 
energy levels. 

(iii) At low temperature, we consider polaron formed 
by single impurity molecule immersed in the Fermi sea. 
We find that the polaron binding energy is larger com- 
pared to a wide MFR at the BEC side and the resonance 
regime, while it is smaller at the BCS side. 
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Two-body Problem. Using D as length unit r — > r/D 
and Eu = h 2 /(mD 2 ) as energy unit E — > E/Eo : the 
two-body Schrodinger equation can be written in dimen- 
sionless form 



1. 



1 -3cos 2 6> 



<5 = E^ 



(2) 



in the regime r > r c and the boundary condition is 
given by ty(r = r c ) = 0. Expanding the wave func- 
tion in terms of spherical harmonics as \&(r, 9, tp) — 
J2i m lRim(r)Y lm (9,(p), the Schrodinger equation can be 
reduced to 



2dr 2 



Rlr, 



1 



2 H%,(r)R l , m {r) = ER. 



Irn 



(3) 



where Jf<™> = l{l+l)6 u , + (Y lm \(l-3cos 2 d)\Y l>m ) /r. H w 
couples all? to I ± 2 and m is still a good quantum num- 
ber. Moreover, since the dipolar coupling decays (1/f 3 ) 
faster than the centrifugal barrier term (1/r 2 ), at large 
distance, coupling between different I effectively vanishes. 
The asymptotic behavior of the scattering wave function 
still can be decoupled as different partial wave channels, 
so one can introduce a scattering length for each partial 
wave channel [TT| . 

Here we will focus on the s-wave channel since it has no 
centrifugal barrier and the interaction effect is expected 
to be the strongest. However, it is easy to note that 
Hqq = 0, which means one cannot obtain any interac- 
tion effect if directly projecting the Schrodinger equation 
Eq. ([2| into the \l,m) = 1 0, 0) state. The interaction 
effect in the s-wave channel comes from its coupling to 
higher partial wave channel. Following a simple second- 
order perturbation argument, one can obtain an effective 
potential from a virtual process between |0,0) and |2,0), 



(00|(1 -3cos 3 0)|2O)(2O|(1 - 3cos 2 0)|OO)2 



D 2 

■ r 4 • 



Hence, as D increases, the effective potential becomes 
deeper and deeper, which introduces a series of bound 
states and causes scattering resonances. This is the basic 
mechanism of DHR. Since the energies of those inter- 
mediate states for this induced interaction are h 2 l(l + 
l)/(m(r) 2 ) > 2E F for I = 2, 4, 6, . . . , we expect that the 
presence of a Fermi sea and the Pauli blocking effect will 
not strongly affect the induced interaction, which justifies 
using the induce interaction in many-body studies. 

Following Ref. [10], we first diagonalize the matrix Hw 
by a unitary matrix X(r) as X> (r)H(r)X(r) = A(r). 
Translating the radial wave function as <I>; = Xt,Ri', Eq. 
([3]) will become 



1 d 2 , , , d , , . , 
2 + qiv {r)—®i> + v w (r)$// 



2dr 2 



dr 



E§ 



(4) 



where q(r) — ~^X^(r)dX(r)/dr and v(r) = A(r)/r 2 — 
^X'(r)d 2 X(r)/dr 2 . At a large distance, because H(r) 




FIG. 1: (a) Scattering length a s /r c as a function of D/r c ; (b) 
kD cot 6k as a function of (kD) 2 at different D/r c as marked 
in (a), (c) ro/D as a function of —D/a s across different reso- 
nances. 



becomes effectively diagonal, X becomes an identity ma- 
trix and still has the same asymptotic behavior as 
Ri. At an intermediate and short distance, the diagonal 
part of v(r) contains the centrifugal potential contribu- 
tion for all / ^ 0, which separates s-wave from other par- 
tial waves. Because of this energy separation, we can take 
the approximation that only keeps the diagonal terms of 
qui and vw . The infinite number of coupled differential 
equation of Eq. Q is reduced to a set of decoupled equa- 
tions 



E*i 



(5) 



In Eq. ([5]), vqq(t) 7^ is the induced interaction in s- 
wave channel. Numerically solving Eq. ^ for I = 0, we 
obtain its asymptotic behavior of $z = o as sin(fcr + 6k), 
where k = \j2mE jh. 

Scattering Resonances and Effective Range. Once we 
obtain s-wave phase shift from solving Eq. |5]), we can de- 
fine s-wave scattering length as a s — ~ lim^o t&nSk/k. 
We plot a s /r c as a function of D/r c in Fig. 0(a). It dis- 
plays a series of resonances as D/r c increases, and the res- 
onance window in D/r c axes gets wider for larger D/r c . 
The locations of resonances can be determined by an in- 
tuitive WKB estimation [TU]. This approach precisely 
captures all s-wave resonances with relatively wide win- 
dows [9], and misses some extremely narrow ones which 
are due to a bound state in a higher even partial wave 
channel and coupled to s-wave channel by the ignored 
off-diagonal coupling in qu< and vw . We ignore them be- 
cause they are too narrow to be studied experimentally. 

We expand the phase shift as k cot 5k = — l/a s + 

rofc 2 /2 + In Fig. []Jb) we find that kD cot 6k is a 

very good linear function of (kD) 2 , and the positive slope 
means that tq is positive. In Fig. |T](c), we plot ro/D as 
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a function of D/a s for different resonances with different 
D/r c . It is very remarkable that all these curves across 
different resonances perfectly coincide with each other, 
which shows that r /D is a universal function of D/a s . 
In another word, the short-range physics (r c /D) is irrele- 
vant in determining tq. Inspired by the relation between 
effective range and physical range of a square well poten- 
tial, we fit the curve by tq/D = ap+aiD / 'a s +a2(D / 'a s ) 2 , 
and find a = 1.8390± 0.0006, and a x = -I0.168±0.002 
and «2 = 158.2 ± 0.1. All the high order coefficients 
are several orders of magnitudes smaller. At resonance, 
ro ~ 1.84D. Noting that D can be of the same order of 
1/&F, we can therefore realize resonances in cold atom 
system with a large and positive effective range. This is 
very important of using cold atom to simulate nuclei and 
neutron matter [T31 HI] ■ 

High-Temperture Regime: In the rest part of the paper, 
we will investigate the many-body effects of a positive 
effective range. We shall compare DIIR with a zero-range 
model for a wide MFR, which gives a constant phase 
shift kcotSk = — l/ct s for all k, and ro = 0. First, we 
study the high-temperature regime by second order virial 
expansion, because resonant interaction manifests itself 
in interaction energy even above degenerate temperature, 
and can be measured by spectroscopy method easily. 

By second order virial expansion, we have 



\E b \/(k h T) 



dk dSk 



(6) 



/ tt dk 

where we have ignored all the contributions from high 
partial waves, and A = ^2irh 2 / (mk^T). And the inter- 
action energy is given by 



Eint 



h_ V2 T db2 
V2 3 dT 



(7) 



In Fig. [2]ja), we compare tint Akin for a DIIR with a wide 
MFR, and in Fig. |2fb), we show e; n t/ekin for a DIIR with 
different temperatures, where ekin = ik^Tn(nX 3 )/2. In 
the side with a s > 0, negative e; n t includes the contribu- 
tion from bound state, while positive ei n t excludes bound 
state contribution, which corresponds to the physics of 
" the upper branch" . 

First, for scattering states, it is known that, ei n t/ekin 
approaches =f0.5 as one approaches a wide MFR from 
negative and positive side of a s [TT] . While for a DIIR, 
we find for a s < 0, the attractive interaction energy is 
weaker for a DIIR, and for a s > 0, the repulsive interac- 
tion energy is stronger. This is strongly in contrast to a 
narrow MFR (r < 0) where the opposite effect is found 
by Ref. [16j . This can be understood from an argument 
based on energy dependence on scattering length. Let us 
consider a situation that a s not very close to resonance, 
and ro is non-zero, but is still small compared to A 2 /a s . 
Because a(k) — l/(l/a s — r k 2 /2), (k denotes a thermal 
average of k, and is of the order of 1/A), for a s < 0, 
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FIG. 2: (Color online) (a) eint/ekin as a function of —X/a s for 
a DIIR (solid line) and a wide MFR (dotted line). X/D = 1. 
Inset: same plot for \/D — 70. (b) eint/ekin for DIIR with 
different temperatures. X/D — 2 for green blue dashed line; 
X/D — 1 for blue solid line and X/D — 0.5 for red dotted 
dashed line. Positive eint/ekin at positive a s side excludes the 
contribution from bound state. 



one finds \a(k)\ < \a s \ for ro > 0, and \a(k)\ > \a s \ for 
r < 0; And for a s > 0, one finds \a(k)\ > \a s \ for r > 0, 
and |a(fc)| < \a s \ for r < 0. This is consistent with ro 
dependence of interaction energy shown in Fig. [2ja). 

An alternative way to understand this result is through 
analyzing energy levels. Let us consider two-atom con- 
fined in a large hard sphere so that the two-body wave 
function ip(r) has to vanish at r = R. The energy lev- 
els of s-wave states is determined by E n = f^k^/m, 
where k n satisfies the equation k n R + dk n = nn. In 
the limit a s — > 0~, Sk = and fc° = nir/R; while for 
a s — > + , 8k = tt and fc" + = (n — \)tt/R. For a wide 
MFR, at resonance, 8k = 7r/2, and all k n is given by 
(n — l/2)n/R, which is responsible for e^t/ckm = ±0.5 
at resonance [16] . While for a DIIR, at resonance, k satis- 
fies kR + arctan(2/(rofc)) = nix. At a given temperature, 
the system is dominated by states with k ~ 1/A. Hence, 
arctan(2/(rofc)) is always smaller than 7r/2 if ro > and 
it decreases to zero as temperature increases. Therefore, 
the downshifts of the energy levels E® — E n are smaller 
compared to a wide MFR, and they decrease as temper- 
ature increases. And the upshifts of the energy levels 
E n — E® are larger compared to a wide MFR, and they 
increase as temperature increases. This is also consis- 
tent with temperature dependence of ei n t/ekin as shown 
in Fig. gb). 

The effective range r also has effect on the bound state 
energy. The bound state energy is given by the pole of 
scattering amplitude f{E) — l/(l/a s — r E/2 — i^/E), 
which gives 



E 



h 2 1 - i/l - 2r /c 
mo 2 r /a s 



ma- 



2a s 



(8) 



For the same a s , the binding energy l^l is larger with 
positive ro. For same reason, for a narrow MFR with 
negative r , the binding energy is smaller. In the inset 
of Fig. |2^a) we show that at a lower temperature, the 
attractive interaction energy of a DIIR will finally ex- 
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FIG. 3: (Color online) (a) Polaron energy £/Ef as a function 
of — l/(fcpa s ) for kpro = (wide resonance, solid line); kpro = 
0.2 (dashed line) and k?ro — 0.5 (dotted line); (b) S/Ep as a 
function of kpro for l/(fcpa s ) = —1 (solid line); l/(fcpa s ) = 
(dashed line) and l/(fcpa s ) = 0.2 (dotted line) 



ceed that of a wide MFR. This is because for a lower 
temperature, the contribution of the lowest bound state 
will gradually become dominative, and the larger binding 
energy finally overwhelms the smaller contribution from 
the scattering states. 

Single Impurity Problem at Low Temperature: In gen- 
eral, the low-temperature many-body physics at res- 
onance is too sophisticated to be studied by simple 
method. Here, in order to illustrate that the effect of 
ro at low temperature could be different from that in 
high temperature, we consider a simpler situation where 
one single impurity molecule is dressed by a Fermi sea of 
majority molecules and form a polaron. For wide reso- 
nance, such a problem can be well described by a vari- 
ational wave function or T-matrix approach |18j . Here 
we apply similar approach to DIIR and obtain following 
self-consistent equation for polaron energy £ 



|q|<|k F | 



1 



Yin. 

47T 



I{£) 



(9) 



where J(£,q) = £ |k|>fcF k -e^s 

\ k 1 p+q— k P 

E| k |< fcF i> and el = k 2 /m, e£ = k 2 /(4m). The dc- 
tails for obtaining this equation is similar as the same 
impurity problem across a narrow MFR 19] . We find 
that the effective range has more dramatic effect on po- 
laron binding energy at the BEC side than at the BCS 
side, as shown in Fig. |3][a). And we find that at the 
BCS side, a positive ro will decrease the polaron bind- 
ing energy \£\, while it will increase |<£" | at the resonance 
regime and the BEC side. Such an effect is also opposite 
to polaron nearby a narrow MFR [19j . and for narrow 
resonance, such effect has been observed in a recent ex- 
periment [20]. In both cases, it is caused by and can 
be understood from the energy dependence on scatter- 
ing length [19]. It will be very interesting to observe the 
counter effect in dipolar gases of polar molecules, so that 
the effects of effective range in many-body systems will 
be established in a comprehensive way. 
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